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,  ABSTRACT 

H 

The  possible  scaling  of  the  axisymmetric  component  of  spectral 
density  of  wall  pressure  in  turbulent  flow  along  a  slender  cylinder 
is  reexamined  with  use  of  a  formal  low-wavenumber  expansion  in 
terms  of  fluctuating  Reynolds  stresses  as  sources.  A  rudimentary 
model  of  the  sources  is  formed  that  conforms  to  the  principle 
of  local  similarity  near  the  wall  and  to  probable  implications  of 
fluctuating  and  mean  velocity  profiles  measured  in  cylindrical 
boundary  layers.  To  the  extent  that  wall  pressure  is  exclusively 
due  to  such  a  turbulence  field,  to  zero  order  in  wavenumber,  at 
least,  it  is  indicated  to  scale  with  cylinder  radius  (a)  and  to 
agree  in  functional  dependence  with  the  current  model  for  this 
low-wavenumber  component.  If  the  turbulence  field  (despite  lack 
of  experimental  evidence)  involves  in  addition  a  weak  component 
that  scales  with  boundary-layer  thickness,  then  at  sufficiently 
low  frequencies  and  wavenumbers  the  corresponding  contribution 
to  wall  pressure  (when ($/a>>l)  may  become  dominant,  and  the  ex¬ 
trapolation  of  the  currentfnodel  to  very  small !j^/U<)o_and  the 
associated  prediction  of /diameter  dependence  may/ then  fail,  "''n 
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1.  INTRODUCTION 

The  question  of  accurate  and  validated  modeling  of  the  axisym- 
metric  component  of  turbulent  wall  pressure  on  a  towed  flexible 
cylinder  remains  one  of  intense  interest.  Apart  from  elements  of 
uncertainty  in  the  parameter  domain  of  common  concern,  it  is 
becoming  important  to  extrapolate  into  diameter/speed/frequency 
domains  where  the  experimental  basis  for  modeling  is  more  indirect 
and  less  secure.  Such  modeling  is  needed  to  predict  prospective 
perfoxmance  in  future  innovative  as  well  as  more  conventional 
improved  systems  and  thus  to  guide  development. 

In  this  regard  the  question  particularly  arises  as  to  whether 
a  current  tentative  model  for  the  wavenumber-frequency  spectral 
density  of  the  subject  axisymmetric  harmonic  of  wall  pressure, 
Pg(k,u),  incorporating  as  the  sole  pertinent  length  scale  the 
cylinder  radius,  a,  is  likely  to  provide  a  correctly  indicative 
result  when  extrapolated  to  very  small  values  of  wa/U^.  This 
memo  attempts  further  discussion  of  this  topic  by  reference 
especially  to  available  measurements  and  further  plausible  model¬ 
ing  of  turbulent  flow  at  large  values  of  6/a,  the  ratio  of 
boundary-layer  thickness  to  cylinder  radius.  This  represents 
a  line  of  development  based  in  part  also  on  previous  formal 
and  8emiempirical  modeling  pursued  in  the  ONR  program  (see 
especially  Refs.  1-3). 

A  more  complete  application  of  current  modeling  to  measured 
noise  in  high-speed  towing  tests  of  quiet  modules  than  that  under¬ 
taken  in  Ref.  2  is  also  needed.  In  particular,  the  preliminary 
analysis  of  data  from  the  LSM  used  a  simplified  and  somewhat  in¬ 
adequate  representation  of  the  wavenumber  filtering  action  of  the 
specially  mounted  hydrophone  employed.  Such  further  application  of 
modeling  may  be  undertaken  in  towed-array  exploratory  development 
outside  the  scope  of  the  more  basic  ONR  program  and  would  be  com¬ 
plementary  to  the  present  work.  An  updated  appraisal  of  the  residual 
noise  level  attainable  by  towed  arrays  could  then  also  be  performed." 
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•  A  related  question  not  addressed  further  here  is  whether, 
for  given  speed  and  streamwise  position  along  a  cylinder,  there 
is  a  diameter  below  which  the  flow  is  no  longer  fully  turbulent. 
Suggested  criteria  for  such  a  diameter  (or  related  Reynolds 
number)  have  been  reviewed  in  Ref.  30  (p.  124),  but  later  mea¬ 
surements  discussed  in  Ref.  19  (pp.  22-23)  indicate  that  a 
transition  diameter,  if  it  exists  at  all,  may  be  exceedingly 
small.*  Below  such  a  transition  threshold,  the  present  semi- 
empirical  modeling  and  resulting  scaling  of  wall  pressure  for 
turbulent  flow  would  not  apply. 


2.  ELEMENTS  OF  PREVIOUS  MODELING  AND  THE  POSSIBLE  ROLE  OF 
DISTINCT  LENGTH  SCALES 

Initially,  attention  is  restricted  to  the  form  of  Pg(k,u>) 

at  negligible  wavenumber,  i.e.  P  (0,u>).  According  to  the  current 

0 

model 

P  (0,w)  ■  Cp2v*a“3  (ui/v*+h/ba)“5,  (1) 

o 

where  C  and  h/b  are  numerical  constants  and  v*  is  identified  as 
the  wall  friction  velocity  and  taken  proportional  to  U^,  say 


*For  a  cylinder  in  water  at  a  flow  speed  of  15  kt,  for  example, 
it  is  indicated  that  the  transition  diameter  is  $  0.01  cm. 
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v*  =  0.041^.  Though  this  form  was  suggested  by  rudimentary 

modeling  considerations/  it  is  supported  mainly  by  inference 

from  certain  noise  measurements  within  liquid-filled  towed-array 

modules.  Indications  have  been  adduced  that  the  appropriate 

length  scale  in  P  Uc,u>)  is  cylinder  radius,  as  implied  by  Eq.  1, 

o  , , 

and  the  preliminary  analysis  of  measurements  in  the  LSM,  a 

module  of  larger  radius,  indicated  rough  consistency  with  the 

notion  of  a  as  the  length  scale  and  with  the  model  form  of 

P  (k ,u>)  that  underlies  Eq.  1  (see  Sec.  3).  Nevertheless, 
o 

these  indications  are  not  yet  uniquely  compelling,  and  the  form 
(1)  presently  rests  mainly  on  a  single  module,  XN-55.24 

Model  (1)  or  any  model  of  Po(0,u>)  based  on  a  single  length 
(and  velocity)  scale  has  the  form4 

Po(0,o>)  «  p2v4a2F(o>a/v*)  .  (2) 

From  an  uninformed  point  of  view,  it  may  naturally  be  asked 

whether  in  a  more  exact  model  there  may  not  be  included  additively 

in  P  (0,o>)  a  similar  second  term,  p2v*62G(u>S/v*) ,  say,  that  scales 
o 

rather  on  boundary-layer  thickness.  Suppose,  as  a  conjectural  example, 
that  these  two  functions  F  and  G  each  had  the  dependence  implied 
in  Eq.  1,  so  that  one  might  write 

P  (0,u>)  ■  p2vi[C  a"s(w/v*+h  /b  a)"*+C  6"s(u)/v*+h  /b  S)“5].  (3) 

0  I  11  2  "22 

If  the  respective  intensity  and  scale  coefficients  satisfy  the  con¬ 
ditions 

C  »(a/«)*C  but  C  (b  /h  )5<<(6/a)2C  (b  /h  ) 5 
1  2  111  222 

(as  would  be  so  if,  for  example,  6>>a  and  C  ~C  ,  b  ~b  ,  h  ~h  ) ,  then 

121212 

the  4-scaled  contribution  would  become  dominant  at  sufficiently 
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small  wa/v*  even  though  the  a-scaled  contribution  was  dominant 
at  greater  values,  including  even  a  plateau  interval,  as  shown 
in  the  sketch  below.  In  such  event,  extrapolation  to  lower 


diameter-frequency  products  on  the  basis  of  a  model  fit  to  mea¬ 
surements  in  the  range  of  dominance  by  the  a-scaled  contribution 
could  grossly  underestimate  the  noise  below  the  steep  rise  due 
to  the  6-scaled  contribution.  In  particular,  the  dependence  of 

P  (0,u>)  on  a  as  a2  for  b  wa/h  v*<l  predicted  from  the  a-scaled 
0  1  1  ~ 

contribution  alone  would  be  replaced  by  independence  of  a  at 
values  of  a  small  enough  that  the  <S-scaled  part  became  dominant. 

The  preceding  hypothetical  account,  though  illustrative, 
is  oversimplified  and,  if  valid  at  all,  requires  modification. 

In  particular,  even  a  contribution  to  axisymmetric  wall  pressure 
associated  with  a  fluctuating  velocity  that  scales  with  6  may 
involve  also  cylinder  radius,  on  account  of  its  geometric  role, 
contrary  to  the  form  of  the  second  term  in  (3).  Hence,  in 
addition  to  delineating  and  modeling  possible  scaling  properties 
of  the  turbulent  boundary  layer  in  the  light  of  present  under¬ 
standing  of  wall  turbulence  and  measurements  in  cylindrical 
geometry,  it  is  required  to  invoke  an  explicit  formal  relation 
of  axisymmetric  wall  pressure  to  the  turbulent  velocity  field. 
This  is  pursued  in  Sec.  4. 
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If  a  model  spectrum  with  additive  terms  characterized 
respectively  by  turbulence  scales  a  and  6  (with  additional 
geometric  dependence  of  both  terms  on  a)  is  a  viable  possibility/ 
there  is  no  suggestion  of  it,  at  least,  in  the  profile  of  the 
rms  stream  wise  component  u  of  fluctuating  velocity  in  the 
boundary  layer  on  a  thin  cylinder  measured  by  Leehey  and 
Stellinger.s  (See  Fig.  1,  repeating  Ref.  6,  Fig.  2.3-6,  based 
on  Ref.  5,  Fig.  8.)  This  profile,  apart  from  the  rise  in  the 
viscous  sublayer  and  a  peak  somewhat  outside  it,  conforms  well,  as 

6  — k 

shown  and  analyzed  by  Kronauer,  to  dependence  as  r  ,  i.e. 
w 

u«v*(a/r)  ,  all  the  way  to  r=26a,  while  in  the  subject  case 
the  edge  of  the  boundary  layer  falls  at  r-44a.  Clearly  (though 
even  prior  to  seeing  the  measured  profile  it  would  perhaps  be 
naive  to  imagine  the  turbulent  dynamics  to  permit  it)  there  is 
not  a  contribution  to  the  profile  («r  that  extends  only  to 
a  distance  that  is  a  modest  multiple  of  a  and  is  associated  with 
turbulent  motion  with  scale  ~a,  followed  at  larger  wall  distances 
by  a  distinct  contribution  that  extends  to  r~<S  and  is  associated 
with  a  turbulent  motion  of  scale  ~S, 

On  the  other  hand,  despite  the  cited  monolithic  profile,  a 
contention  might  be  made,  along  the  lines  of  the  discussion  by 
Bradshaw,7  following  Townsend,8  referring  in  their  case 
to  planar  turbulent  boundary  layers  with  or  without  mean  pressure 
gradient,  that  the  turbulence  is  nevertheless  composed  of  two  such 
components  of  disparate  scales.  The  large-scale  wake-like  com¬ 
ponent,  described  as  "inactive,"  would  be  thought  not  to  transfer 
much  momentum  nor  to  interact  appreciably  with  the  "universal" 
component,  which  would  scale  rather  with  distance  y  from  the  wall 
when  y<a.  Part  of  the  modeling  introduced  below  is  intended  to 
reflect  this  possible  contention. 
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3.  SELECTIVE  REVIEW  OP  MODELING  PRESSURE  ON  A  CYLINDICAL  WALL 

Elements  of  previous  modeling  will  now  be  further  outlined. 

In  Ref.  9  (but  based  on  the  earlier  Ref.  10)  a  type  of  model 

was  advanced  for  the  azimuthal  harmonic  spectral  densities  P  (k,w) 

m 

(m  integer)  of  turbulent  wall  pressure  on  a  cylinder  by  writing 
t  (m+%)  /a 

Pm(k,uj)  =  I  dkaP(k,ks,o>)  (4) 

(m-)5)/a 

and  choosing  a  model  spectral  density  P(k,k3,w),  which  depends  on 
a  second  wavenumber  analogous  to  the  transverse  component  in 
the  instance  of  a  planar  boundary.  For  the  planar  case  a  theoreti¬ 
cally  attractive  model  for  the  wavevector-frequency  spectral 
density  of  wall  pressure  in  the  low-wavenumber  domain  that  con¬ 
formed  as  well  as  any  to  the  laboratory  measurement  of  low- 
wavenumber  planar  wall  pressure  by  Jameson12  was  the  scale- 
independent,  wavenumber -quadratic  one  of  the  form 

P(k  ,k  ,u>)  »  C  p2viu>-3  (v*K/w) 2 ,  K2  =  k2+k2  (5) 

13  2  13 

(C2  constant).  Eqs.  4  and  5  yielded  for  application  in  some 
domain  restricted  at  least  by  the  condition  ksui/u^, 

P#(k,w)  »  C2P2vi<i)"5a" 1  (k2+^j-  a"2)  .  (6) 

A  trial  model  to  encompass  the  full  wavenumber/frequency 
domain  of  interest,  including  the  convective  one,  was  constructed, 
in  effect,  as  follows.  The  spectral  density  of  wall  pressure  for 
a  frame  convecting  downstream  at  some  speed  uc  is  obtained  from 
(6)  by  replacing  u  by  u-u  k.  Suppose  that  in  that  frame,  for 
suitably  chosen  u  ,  the  decorrelation  of  axisymmetric  wall 

C 
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pressure  due  to  streamwise  separation  £  is  comparable  with  that 
due  instead  to  time  delay  t  if  ;  and  t  are  related  by  a  dis¬ 
persion  velocity  scale  ~v*,  i.e.  if  5~v*t,  and  suppose  that,  for 
separation  and  time  delay  both  nonvanishing,  the  decorrelating 
effects  add  in  quadrature  (see  also  Refs.  11,  21).  If  the  de- 
correlation  is  also  scale-independent,  we  have  then  to  replace 
w/v*  in  Eq.  6  by 

k+  =  [  (u)-uck)2/(hv*)2+k2],s  ,  (7) 

where h  is  of  the  order  of  unity.  For  sufficiently  low  frequencies 
|u)-uck|  and  wavenumbers  k,  however,  a  cutoff  of  the  rise  in  the 
pressure  spectrum  must  be  effected  by  the  finite  scale,  say  A, 
of  the  largest  contributing  turbulent  motions.  Originally,  it 
was  assumed  for  convenience  that  this  cutoff  wavenumber.  A**1, 
could  be  added  in  quadrature  with  the  frequency-augmented  wavenumber 
(7) ,  so  that  the  generalized  model  based  on  Eq.  6  became 

i  -5/2 

P  (k,u>)  =C  p2v|a-1  (k2  +  J*-  a”2)  [  (u>-u  k)  2/(hv*)  2+k2+A-z] 

0  2  -L*  c 

(8) 

When  this  model  was  applied  to  infer  turbulence-induced  noise 
within  a  liquid-filled  hose  and  results  compared  to  various  test 
measurements  most  suitable  among  those  then  available,  the  in¬ 
dicated  magnitude  of  the  level  coefficient  C  in  (8)  proved — 

2 

possibly  fortuitously  but  certainly  encouragingly— to  have  approxi¬ 
mately  the  value  inferred  for  C  in  the  underlying  planar  model  (5) 

2 

on  the  basis  of  Jameson's  measurements.  When  measurements  became 
available  in  quieter  small -diameter  liquid-filled  modules,  not¬ 
ably  XN-55,  especially  at  higher  speeds,  however,  it  was  found 

that  the  inferred  form  of  P  (0,<d)  would  be  appreciably  better 

0 

matched  if  the  model  (8)  was  altered  simply  by  assuming  the  wave- 
number  scale  A-1  to  enter  the  wavenumber  denominator  linearly 
instead  of  in  quadrature,  yielding 
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PQ  (k,w)  =  C2P2v*a_1  (k2+^j  a-2){[  (oj-uck)2/(hvA)2+k2])5+A-1 } 

(9) 

This  form,  though  mathematically  less  convenient  than  (8) ,  seems 
more  likely  also  on  the  basis  of  possible  explicit  source 
modeling. 

The  relation  (4),  if  P(k,k3,uj)  is  identified  with  the  true 
spectral  density  P(k^,k  ,u)  for  a  planar  boundary,  clearly  ensures 
consistency  of  the  Pm(k,u))  with  the  result  required  in  the 
planar  limit  defined  by  S/a-*0,  ka-M°,  oia/v*-**0  ,  in  which  the 
discrete  circumferential  wavenumber  components  m/a  become 
closely  spaced  compared  with  all  pertinent  wavenumbers.  The 
minimum  necessary  limiting  condition  of  correspondence  between 
the  Pm { k , u» )  and  P(k  ,k3,u»)  does  not  seem  obvious,  however,  and 
may  be  much  weaker  than  that  implied  by  Eq.  4.  For  example,  as 
noted  more  explicitly  in  Ref.  2  (footnote,  pp.  15-16)  there  may 
be  no  unique  correspondence  referring  to  a  single,  finite  value 
of  m.  The  original  postulated  relation  (4)  thus  may  not  be 
valid  even  as  a  limiting  approximation,  at  least  when  m  =  0. 
Correspondingly,  the  modeling  analysis  of  Ref.  1  (even  without  the 
specialized  modeling  assumptions  therein)  indicated  no  basis  for 
the  specific  form  a-2 [1/12+ (ka) 2 ]  of  the  parenthetical  factor 
that  entered  Eqs.  6,  8,  and  9.  It  suggested  that  in  trial  modeling 
the  quantity  1/12+ (ka) 2  be  replaced  by  a  somewhat  more  general 
polynomial  in  ka. 

For  reference,  the  trial  model  form  we  presently  employ 
is  given  by 

P  (k,ui)  *  C  p2v*a-1  (a~2--s  k2)  {[  (u>-u  k)  2/(hv*)  2+k2  l^+A-1 } 

0  2  2  c 

with  A  =  ba.* 

♦Apart  from  any  more  general  modification  that  may  prove  necessary, 
presence  of  a  term  linear  in  k  in  addition  to  a-2  +  s2k2  is  not 
excluded. 
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Regarding  the  length  scale  A  entering  P  (k,w) ,  apart  from 

0 

the  as  yet  imperfect  analysis  of  the  implications  of  parallel 
comparisons  with  noise  measured  in  modules  of  different  diameters, 
alluded  to  above,  other  indications  and  suggestions  were  drawn. 
Measurements  of  frequency  spectra  of  wall  pressure  by  small  (not 
axisynunetric)  flush  hydrophones  at  various  positions  along  a 
towed  cylinder  were  made  by  Markowitz.14  As  noted  in  Ref.  14  (see 
also  Ref.  9),  these  conformed  more  nearly  to  scaling  with  radius 
a  than  with  layer  thickness  6,*  and  measurements  by  other  workers 
were  consistent  with  such  a  result  when  6/a>>l.  It  was  proposed 
in  Ref.  16  (see  also  Ref.  9}  that  a  model  wall  pressure  with  a 
single  scale  A,  e.g.  as  in  Eq.  9,  would  suffice,  but  this  scale 
would  depend  on  both  a  and  6,  becoming  *6  in  the  limit  where 
<5/a-*-0  and  increasing  more  strongly  with  a  than  with  6  where 
6/a>>l. 

4.  FORMULATION  OF  MODEL  SOURCE  SPECTRA  FOR  APPLICATION  TO 
WALL  PRESSURE 

Recent  work  in  the  present  program  provides  the  formal  relation 

of  the  subject  pressure  spectrum  P  (k,w)  to  source  spectra  solely 

0 

of  fluctuating  velocity  products  (Reynolds  stresses). 3  This  seems 
likely  to  be  a  safer  approach  (i.e.,  less  likely  to  mislead)  than 
that  based  on  the  Poisson  equation  for  pressure  (in  incompressible 
flow)  .  The  latter  involves  also,  apart  from  spectra  of  two- 
component  velocity  products,  in  connection  with  a  mean-flow  interaction 
term,  a  spectrum  of  (single-component)  fluctuating  velocity  normal 


*The  measured  frequency-spectral  dependences'  were  not  closely  con¬ 
sistent,  however,  with  the  present  model  or  with  general  theoretical 
expectations  for  turbulent  wall  pressure. 
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to  the  wall.  This  approach  was  formulated  in  cylindrical  geometry 
in  Ref.  1  and  pursued  for  a  specific  source  model. 

The  present  approach  via  pure  velocity-product  sources  permits 
straightforward  approximation  of  the  pertinent  Green's  function 
required  to  exhibit  wall  pressure  explicitly  in  terms  of  assumed 
sources  only  in  the  nonconvective  tail  (U^ 1 1c  1  <tu) .  An  explicit 
result  lias  been  given  only  where,  in  addition,  |k)5<  1 .  Resulting 
modeling  may  guide  expectations,  however,  also  outside  this 
domain.  At  the  present  time,  it  should  be  mentioned,  formal 
results  of  the  respective  approaches  in  this  domain  exhibit  apparent 
differences  that  remain  to  be  elucidated. 

The  relations  from  Ref.  3  to  be  used  are  available  in  a  form 
sufficiently  general  to  include  slightly  compressible  flow.  This 
generalization  is  worth  having  though  beside  the  point  of  the 
present  considerations. 

At  k=«0,  i.e.  to  lowest  order  in  k,  axisymmetric  wall-pressure 

A  /V 

amplitude  is  related  to  amplitudes  Tr£ (r,k,w) ,  T00  (r,k,u>)  (with 

k*0)  of  the  self  products  of  radial  and  of  azimuthal  fluctuating 
velocity  at  radius  r  (wall  distance  r-a) ,  according  to  Ref.  3, 

Eq.  62  (or  Ref.  1) ,  by 

00 

Po(0,w)  *  pj  drr"1 (Tr®-TQ®)  .  (10) 

a 

If  we  imagine  the  higher-order  terms  in  k  also  to  be  included,  it 
is  sensible  to  retain  the  argument  k  of  the  T^  rather  than  setting 
k=0.  In  fact,  since  we  can  only  try  to  find  support  or  motivate 
modification  of  current  model  forms— -not  pursue  rigorous  expansions — 
and  we  wish  to  do  this  with  regard  to  their  full  domain  of  hypotheti¬ 
cal  applicability,  we  feel  free  to  model  the  spectra  of  T^  as 
though  encompassing  also  the  convective  domain. 

The  spectral  density  at  k  *  0  is  given  by  Eq.  10  as 

00 

Po(0,w)  ■  p2j  drr”1!  drr 'r'"lEo (r,r '  ,0,w)  ,  (11) 

a  a 

10 
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where  E0  (r,r '  ,k,to)  denotes  the  (cross)  spectral  density  associated 
with  Tr£-Te°„  (Index  zero  here  refers  to  m.)  In  accord  with  the 
preceding  discussion,  we  wish  to  consider  the  right  member  of  (11) 
also  for  general  k  and  denote  the  result  by 

GO  00 

Pj(k,w)  -  p2J  drr-1 J  dr'r,-1Eo  (r,r '  ,k,u>)  .  (12) 

a  a 

First,  however,  we  consider  the  formation  of  a  plausible 
type  of  model  for  spectra  E^1  (r,r '  ,k,w)  of  the  streamwise  fluctuat¬ 
ing  velocity  component.  The  corresponding  profile  of  the  single¬ 
point  rms  value. 


E  (r)  =  E  Em  (r)  ,  where  E®  (r)  =  [du>fdkEm  (r,r,k,w)  , 
i i  _ _ _  ii  ii  j  I  i  i 

m  J  J  (13) 

has  been  measured,5  as  noted  in  Sec.  2.  This  is  intended  to 

guide  analogous  modeling  of  velocity-product  spectra,  notably 

E  (r,r',k,u)  that  appears  in  Eq.  12.  The  model  is  required  to 
o 

claim  validity  even  where  r,r'  >a  when  6/a>>l. 

We  consider  a  model  characterized  by  a  single  length  scale, 

A,  but  spectra  may  also  be  formed  by  "incoherent"  addition  of  such 

spectra  having  different  scales  (cf .  Sec.  2) .  In  the  spirit  of 

our  usual  modeling  assumptions  we  suppose  that  E°^  (and  likewise 

velocity-product  spectra  such  as  E  )  depends  on  u>  and  k  only  via  k,  , 

0  • 

defined  by  Eq.  7.  We  may  then  write 

E°  i  (r,r  '  ,k,oj)  *  v*A2  (a/r)  **  (a/r  ' )  °  i  (k+  (rr  ' )  1,s,k+  (r  ’-r)  ,k+A) . 

(14) 

The  profile  factor  (aVrr')*5  is  required  to  conform  to  the  measured 
dependence  of  E^tr)  as  a/r  (for  r <6  except  near  the  wall,  Fig.  1), 
at  least  for  the  class  of  functions  Fj  (and  their  generalization, 
F®  ,  to  m^O)  considered  here. 
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A  suggestive  approach  to  a  type  of  model  for  (especially 
when  6>>a)  may  be  based  on  imagining  first  a  region  of  more  or 
less  homogeneous  turbulence  of  fixed  scale  A  without  regard  to 
presence  of  a  physical  cylinder.  The  fluctuating  velocity  field 
may  then  be  imagined  expanded  in  waves  of  definite  k  and  u>  and 
angular  harmonics  m  defined  for  a  cylindrical  coordinate  system. 

If  a  model  is  assumed  for  wavevector -frequency  spectra  of  fluctuat¬ 
ing  components  appropriate  to  the  hypothetical  turbulent  field,  this 
may  thus  be  expressed  by  explicit  transformation  in  terms  of 
spectral  components  characterized  by  m,k,u>,r,r ' . 

To  treat  the  actual  case  where  a  hard  cylinder  is  present 
and  responsible  for  generating  the  turbulent  field  by  its  motion, 
the  essential  step  is  to  permit  the  scale  A  in  the  formally  con¬ 
structed  model  to  depend  on  the  distances  y  =  r-a,  y'  =  a ' -a  from 
the  wall.  For  y»y'«a  (but  outside  the  sublayer),  in  particular, 

7,17 

the  usual  similarity  principle  for  a  planar  wall  implies 

that  A«y.  (We  look  apart,  for  the  moment,  from  any  weak  inactive 
component  scaling  otherwise.) 

Proceeding  along  these  lines,  as  shown  somewhat  more  explicitly 
in  Appendix  1,  an  indicative  class  of  models  may  be  characterized 
(leaving  aside  the  choice  of  A)  by  regarding  the  function  F°t  in 
Eq.  14  as  given  roughly  in  opposite  limiting  domains  by 

(1  7  kL(zr')h<l 

(15) 

L  ”1  1* 

B[kA(rr'P]  ,  kA(rr’)*>l  , 

kA  =  <ki+b2A“2)*  ,  A  =  [  (r'-r)  2+b2A2]J*  ,  (15a) 

and  b  ,  b  are  constants.*  The  function  S,  (z)  in  (15)  approaches 
0  1  i 

a  constant  (of  the  order  of  unity)  for  arguments  z<l  and  decreases 
at  some  rapid  rate  for  z*l. 

^Either  ba  or  bj  could  be  absorbed  in  the  definition  of  A,  but  it  is 
convenient  for  consideration  of  models  where  b^O  or  b|  0  to  preserve 
their  explicit  appearance. 


12 


TM  21 


Chase  Inc 


On  the  basis  of  properties  of  a  function  found  in  Ref.  11  to 
fit  measurements  by  Bullock,  Cooper,  and  Abernathy18  of  cross¬ 
spectra  (y'-y^O)  of  streamwise  turbulent  velocity  over  an  effec¬ 
tively  planar  wall,  we  assume  that  b(  ■  0  in  (15a) .  (A  certain 

illustrative  model  with  b  j*0  but  b  *0  is  also  considered  to  some 

o  1 

extent  in  Appendix  1.)  The  requisite  estimates  based  on  the 
model  are  greatly  facilitated  thereby. 

We  turn  to  modeling  of  velocity-product  spectra  such  as  E(  in 
Eq.  12.  In  view  of  the  quadratic  dependence  of  these  products 
on  individual  velocity  components  and  the  expectation  that  their 
spectra  resemble  that  of  streamwise  velocity  in  their  dependence 
in  limiting  domains  of  the  arguments  (e.g.  for  kAAr<l  or  >1),  we 
assume  that  model  equations  exactly  analogous  to  (14)  and  (15)  for 
E*x  apply  to  where  the  square  of  the  previous  profile  factor 
must  be  used: 

Efl(r,r',k,u>)  -  v*A2(a/r)  (a/r')F#  (k+trr’)1*,  k+(r’-r),  k+A)  ,  (16) 

where  F  is  characterized  again  as  in  (15)  with  S  replaced  by, 

0  i 

say,  S  . 

0 

For  present  indicative  estimation  of  limiting  dependences 
we  interpolate  in  r,r*  between  opposite  limiting  domains  of 
k^  (rr ' )  **  in  (15)  by  assuming  a  convehient  product  form 

F#  »  (1+ak^r)  _ls(l+okAr ' )  "’,*So  (k^Ar)  (17) 

la2"^”1).*  More  significantly,  we  wish  to  assume  also  a  factorable 

form  for  S  .  In  this  connection,  we  note,  in  S  referring  to  the 
0  1 

•Possible  different  presumptions  include  that  F»[l+okA (rr')1]  S0. 
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spectrum  of  (linear)  streamwise  velocity,  the  earlier  fit11  to 

the  measurements  of  Ref.  18  corresponded  to  taking  A  *  (yy1)**, 

■  /?  (and  k^-k+)  (in  15a),  so  that  Ar»(y2+y ,2) if  the 

same  assumption  is  appropriate  for  ,  referring  to  a  spectrum 

of  velocity  products,  and  if  one  had  S  (k^A  )«  exp[-(8k^A  )2] 

o  +  r  +  r 

(8-const.),  for  example,  then  S  “  exp[-(8k.y)^]exp[-(8kJ.y') 2],  a 
product  form.  Actually,  the  fitted  did  not  have  such  a 
Gaussian  form.  Still,  it  is  perhaps  indicative  of  limiting 
dependences  to  assume  a  product  form: 

So(k+Ar}  "  S(k+y)S(k+y')  ,  (18) 


where  the  function  S(z)  begins  to  decrease  rapidly  for  z  greater 
than  a  value  of  the  order  of  unity.  It  is  tentatively  assumed 
also  that,  for  all  k  and  to,  the  source  spectrum  is  negligible 
outside  the  boundary-layer  thickness,  so  that  (notwithstanding 
its  single  indicated  argument)  the  function  S(k+y)  =  0  for  y>6. 
Finally,  a  factorable  form  is  assumed  also  for  A2: 


A2  »  L(y)L(y')  , 


(19) 


where  L  has  dimensions  of  length. 

Eqs.  12,  16,  17,  and  18  yield  for  the  "generalized"  zero-order 
contribution  P°  (k,<o)  to  the  wall-pressure  spectrum 


P°  (k,u>)  «  p2v*a2I2  (k,io) 

I  (k,<o)  «  f  dy[l+ak+ (a+y) 


]-Js (a+y)  “2L  (y) S (k+y)  . 


(20) 

(20a) 
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5.  SPECIFICATION  OF  SCALING  OF  SOURCE  SPECTRA  AND  RESULTS 
FOR  MODEL  WALL-PRESSURE  SPECTRA 

In  the  context  of  this  memo  a  central  question  is  how  to 
model  the  scale  A  (or,  under  assumption  (19),  L)  .  Apart  from 
suggesting  a  favored  rough  dependence  for  A,  we  wish  to  explore 
the  consequences  of  alternative  assumptions.  If  attention  were 
confined  to  y,y'<a,  as  noted  above,  for  consistency  with  the 
principle  of  local  similarity  we  should  take  for  A(y,y')  a 
homogeneous  symmetric  function  of  first  degree  in  y,y',  so 
that  A(y,y)«y.  In  modeling  for  a  planar  wall  in  Ref.  11  it  was 
proposed  that  A«(yy')  .  For  y,y'>a,  looking  apart  from  any 
possible  "non-similar"  component  scaling  on  6,  the  notion  of  a  scale 
«y  (when  y'«y)  loses  any  claim  to  sensibility, 7  and  the  reduced 
influence  of  the  wall  as  its  distance  y  becomes  large  compared  with  the 
wall  radius  appears  to  imply  that  A  should  be  taken  a  homogeneous 
function  of  first  degree  in  a,y,y*  that  increases  with  y  (for 
y'«y)  more  slowly  than  y. 

We  pursue  as  a  representative  model  along  this  line  one  where 
L(y)  is  proportional  to  a  "mixing  length."  This  is  conventionally 
defined  such  that  its  product  by  the  local  shear  velocity  at  y 
yields  the  "eddy  viscosity",  e;  this  in  turn  is  conventionally 
defined  such  that  its  product  by  the  local  mean  velocity  derivative 
U' (y)  in  the  cylindrical  flow  (and  by  density)  yields  the  local 
shear  stress.  I.e., 

e  -  L (t/ p) 
t  *  peU'. 

In  the  approximation  of  circumferential  integrated  shear  force 
independent  of  wall  distance,  which,  as  discussed  in  Ref.  6,  is 
supported  by  the  measured  fluctuating  velocity  profile  measured 
in  Ref.  5,  one  has 


(21a) 

(21b) 


t  ■  pv*/(l+y/a) . 
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(22) 
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Combining  Eqs.  17  and  18  yields 

L  «  (l+y/a)~*v*/U*  (y)  .  (23) 

The  profile  U(y)  for  cylindrical  flow  outside  the  near-wall 
region  (in  the  analog  of  the  log  region  for  planar  wall)  has  been 
measured  extensively  (in  certain  domains  of  Reynolds  number)  by 
Willmarth  et  all 9  and  Leehey  and  Stellinger . 5  An  indicative  approxi 
mation  to  the  measured  U(y)  is  provided  by  "Rao's  hypothesis" 

(see  present  Appendix  2) ;  in  the  log  layer  this  implies 

Av*  Av#/a 

u’(y)  a  H5F  *  TI+FITraT+y 7aT  (A“2*5)  •  (24) 

In  this  approximation  Eq.  23  becomes 

L  «  0.4ain(l+y/a)  (lty/a)**  (25) 


or,  in  opposite  limits, 


L  -*• 


J  0-4y, 

1  0.4  (ay) ^Und+y/a) , 


y/a«i 

y/a>>l. 


The  latter  form  indicates  that,  apart  from  the  weakly  varying  log 
factor,  the  local  scale  varies  as  the  geometric  mean  of  wall 
distance  and  cylinder  radius.  An  empirical  correction  factor 
suggested  for  the  Rao  profile  in  Appendix  2  would  introduce  in 
the  scale  (25)  a  factor  varying  slightly  with  6/a  but  approximately 
independent  of  y  except  near  the  wall.  The  square  of  this  factor 
would  then  be  adjoined  to  Eqs.  27,  36,  and  39  below. 


TM  21 


Chase  Inc. 


Insertion  of  the  mixing  length  L  of  (25)  into  (20a)  and  use 
of  the  properties  of  S(z)  (see  Eq.  18)  yields  in  opposite 
limiting  domains*  where  for  the  case  of  interest  it  is  assumed 
also  that  a/6<<l: 

(const.  *  k+a<<l 

(26) 

(k+a)“ 5/2  ,  k+a>>l  . 

The  result  of  (20)  then  accords  with  the  limiting  forms  given  by 

the  term  corresponding  to  P°  (k,u>)  in  the  current  a-scaled  model 

0 

(9a)  (i.e.  term  a-2,  only,  in  the  factor  a"2+s2k2).  Otherwise 
stated,  assumption  of  a  particular  interpolation  form 
I~(b"1+k+a)"5/2  consistent  with  the  limiting  forms  (26)  yields 
by  (20) 

P°(k,w)  =C  p2v|a“3[k.+(ba)_1 ]  .  (27) 

It  is  well  to  note  results  for  P°(k,u)  under  alternative 

0 

scaling  assumptions .  if  in  place  of  (25)  L  were  taken  to  vary  as  y 
for  y<a  but  to  become  «a  for  y>a,  the  limiting  forms  (26)  and 
hence  a-scaling  as  in  the  current  model  would  again  apply.  If 
L  were  taken  (quite  unjustifiably)  to  vary  as  y  for  all  y<6, 
the  constant  I  of  (26)  for  k+a<<l  would  become  instead  '-An  (1+6/a) 
in  the  subrange  k+6<<l  and  ~An(l+l/k+a)  for  k+6>>l,  with  corre¬ 
sponding  modifications  implied  in  (27) . 

Suppose,  on  the  other  hand,  that  there  is  an  inactive  con¬ 
tribution  to  fluctuating  velocity,  as  mentioned  earlier,  that 
scales  entirely  with  6  and  that  no  profile  factor  (a/r) (a/r') 

applies  to  this  component  of  E  .  Then  A«6  and,  in  place  of  (18), 

0 

S  ■  S  (k. 6) •  Then  in  place  of  (20)  and  (20a)  we  obtain 
oo  + 

P°  (k,u>)  -  p2v*62I2  (k,u>)Sn  (V  6)  (28) 

0  "0  0  x 

I5(k,u)  ■  /06dy[l+ak+(a+y)  ]“**  (a+y)  ~l .  (28a) 
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In  limiting  domains 


I 


6 


1 


in (1+6/a) , 
in(l+l/ak+a)  , 

(k+a)_J5  , 


k+a<<l,  k+6«l 
k+a<<l<<k+6 

k+a>>l,  k+6>>l 


An  example  interpolation  form  then  yields  from  (28) 


(29) 


P°  (k,to)  =  p2v*52S  (k.  6)  (k  +l/ain2[l+(c,  k, +c,6-1 )  a”1]).1  , 

(30) 

where  c+  and  are  constants.  Setting  k+-*-u>/hv*  for  k*0  yields 

the  corrected  form,  under  present  assumptions,  of  the  hypothetical 

second  term  in  Eq.  3  for  P  (0 ,«) .  In  the  extreme  low-k,  low-u> 

0 

limit  where  k+6<l,  in  particular,  the  resulting  equation  then 
takes  the  form 


P  (0,uj)  -  p2v2[C'a2  +  C'62in2  (I+6/a)  ],  (31) 

o  1  2 

where  C',C'  are  constants. 

1  2 

Such  a  hypothetical  6-scaled  component  of  the  turbulence 
might  have  somewhat  the  structure  of  the  wake  of  an  axisymmetric 
body,  as  suggested  in  Ref.  16.  a  model  for  the  corresponding 
turbulent  sources  of  pressure  might  be  constructed  somewhat  as 
though  this  turbulence  were  homogeneous  (modifying  the  similar 
course  taken  in  Ref.  16  for  wall  pressure  itself);  this  model 
fora  would  differ  somewhat  from  that  taken  in  Eq.  28.  In  the 
absence  of  motivation  by  specific  experimental  results,  detailed 
development  of  this  kind  seems  unrewarding. 

Reverting  to  the  generalized  similarity  model  that  led  to 
Eq.  27,  we  wish  to  consider  also  terms  of  higher  order  in  k  in 

the  low-k  expansion  of  P  (k ,«) .  Eq.  62  of  Ref.  20  for  the  two- 

o  ^ 

term  expansion  of  the  amplitude  p  (k,w) ,  from  which  the  spectrum 

o 

P  may  be  formed,  may  be  written 
0 
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00 

Po  (k,w)-*-p®  (k,u>)  -2pk|  dr{TggU)-1U  '  (r)  [fcn (r/a) +C ' ]-iTr£} , 


where 

C*  -  4n(a»a/Uao)+iln|(»)/kUoo|+y«,n|l-(a)/ck)2|"1+0.116+ia,  (32a) 

p®  (k,u>)  -  pTdrr-1  (Tr®-T@J)  ,  (32b) 

■*a 

a«ir/2  or  0  for  |k|<u»/c  or  >ai/c,  respectively,  and  it  is  understood 
that  (r,k,w)  [cf.  Eq.  10,  in  which  T^-T^  (r,0,w)  ] . 

P°(k,(i>),  we  recall,  represents  the  spectral  density  formed  from 
(32b)  (see  Eq.  12) . 

A 

The  term  in  (32)  containing  in  its  integrand  a  factor  Trx 
has  the  same  form  as  the  term  (32b)  already  considered  except 
that  ^rp“^eQ  is  replaced  by  i2krTrx.  The  corresponding  con¬ 
tribution  to  P  (k,oi),  say  Prx(k,u>)  (excluding  the  cross-term 

/V  o  o 

between  the  Tr£  term  and  other  terms  in  Eq.  32), under  assumptions 
similar  to  those  leading  to  Eqs.  20,  20a,  may  similarly  be  written 
in  the  form 

Prx(k,o>)  =  p2v*a2I2  (k,w)  ,  (34) 

0  £  i  X 

Irx(k»a))  ~  kf  dy[l+ak+  (a+y)]  (a+y)  -1L(y)  srx(k+y)  ,  (34a) 

where  the  dimensionless  model  source  function  Srx  is  related  to  i2Trx  as 
S(k+y)  was  to  $rr”Teg. 

Assumption  of  the  mixing  length  L  of  (25)  as  L(y)  in  (34a) 
and  of  properties  of  Srx(z)  similar  to  those  assumed  above  for 
S(z)  yields  in  opposite  limiting  domains,  where  also  a/6<<l, 


(6/a) 

Irx  ~  le»  *  <  (k+a) 


k+6<<l 

k+a<<l<<k+6 


(k+a)“5/2,  k+a>>l  . 
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3TX 

The  corresponding  limiting  forms  for  P  given  by  (34)  may  be 

o 

subsumed  by  am  assumed  interpolation  to  obtain 

P^x(k,«)~Cp*v*a-*(lca)2[k++(ba)-I]"''[k++(b6S)-l]“1f  (36) 

where  b^  is  a  constant. 

Finally,  we  consider  the  term  in  (32)  explicitly  proportional 

.  A 

to  TQ°.  The  corresponding  contribution  to  P  (k,ut) ,  say  p”(k,u), 
under  the  usual  assumptions,  with  U'(y)  approximated  by  (24)  and 
L  again  taken  as  in  (25) ,  may  be  written 

P® (k,w)  ■  p2v*a2lQ (k,w) ,  (37) 

« 

I0  (k,  u)  »  (kv*/o»)  Jodya_1[l+ak+(a+y)]  (l+y/a)_3/2[in(l+y/a)4C,]S0  (k^y) , 

(37a) 

where  S0  is  related  to  -2T0q  as  S  was  to  Tr£-T00 .  In  opposite 
limiting  domains,  for  a/5<<l, 

II,  k+a>>l 

(38) 

(k+-a)"*/2  ,r  k+a>>l . 

ttf  Eq.  32a  yields  J  C  * | <1 »  in  the  upper  form  here  C'  should  be 
replaced  by  a  coefficient  -l.)  Recalling  that  we  are  dealing  with 
a  formal  expansion  for  MDgB/u<<l  (and  tki6<<l),  implying  a  domain 
where  «ahv*k+,  we  replace  u  accordingly  in  (38),  so  that  the 
result  conceivably  has  some  applicability  in  a  model  encompassing  also 
higher  (e.g.  convective)  values  of  k.  The  resulting  limiting  forms 
given  by  (37)  for  P  may  be  subsumed  by  an  assumed  interpolation 
to  yield 

P®  (k,w)  —  C0 1 C •  1 2p2v*a”  *  (ka)  2k"2[k++(ba)"1  ]”*  ,  (39) 

where  C0  is  a  constant. 
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6.  FURTHER  DISCUSSION  AND  CONCLUSIONS 

To  reiterate,  a  decrease  in  the  spectral  density  of  axi- 
synunetric  wall  pressure  with  decreasing  cylinder  diameter — 
ultimately  as  a2 — is  predicted  at  low  dimensionless  frequencies 
and  wavenumbers  if  the  pressure  in  that  domain  scales  on  diame¬ 
ter  alone  (see  Eq.  2) .  In  the  work  here,  the  likely  most 
plausible  scaling  together  with  rough  modeling  have  been  intro¬ 
duced  to  characterize  the  dominant  component  of  the  turbulent 
field;  at  relatively  small  distances  from  the  wall  the  model 
spectra  conform  to  the  usual  principle  of  local  similarity  and 
are  otherwise  consistent  with  relevant  measurements  for  cylin¬ 
drical  boundary  layers.  It  is  found  that  for  this  assumed 
source  scaling  the  corresponding  wall  pressure  at  negligible 
wave  number,  at  least,  indeed  scales  with  cylinder  diameter 
(see  Eqs.  26,  27) . 

At  the  same  time,  the  predicted  rate  of  decrease  with  diameter 
could  substantially  vanish  in  some  domain  of  very  low  frequency 
and  wavenumber  if  there  is  even  a  weak  distinct  contribution  to 
wall  pressure  due  to  a  component  of  the  turbulence  field  that 
scales  rather  on  boundary-layer  thickness  [e.g.,  see  the  model 
form  (30)  or  its  limit  (31)].  There  appears  no  specific  present 
experimental  evidence  for  such  a  component,  but  an  identifiable 
field  of  this  character  may  presumably  be  present,  though  weak. 

The  relative  contribution  of  this  component  to  wall  pressure 
in  the  low-frequency,  low-wavenumber  domain  naturally  depends  on 
its  relative  source  intensity,  but  also  on  a  factor  somewhat 
like  (6/ a) 2 (e.g.,  see  Eq.  31).  It  may  help  fix  expectations, 
therefore,  to  estimate  6/a  for  given  radius  a  and  effective 
distance  x  from  the  transition  point  upstream.  With  6  as  con¬ 
ventionally  defined  and  v*/U  *  0.04  for  the  flow  domain  of  interest, 
Eq.A2-7  of  Appendix  2,  referring  to  6/a>>l,  gives 

6/a  *  y (x/a) 2/3  , 
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where  the  value  of  y  (which  may  depend  weakly  on  Reynolds  number) 
indicated  by  the  measurements  of  Ref.  19  is  =  0.09  or  by  those 
of  Ref.  5  -  0.13.  For  x/a  =  4000,  e.g.  x  =  50m  and  2a  =  2.5  cm, 
assuming  y  =  0.12  yields  6/a  =  30.  Clearly,  then,  we  will 
typically  be  concerned  with  large  values  of  6/a  such  that 
101og(6/a)2  may  be  -  20. 

With  regard  to  the  opposite  domain  of  high  frequencies  but 
low  wavenumbers,  the  present  plausible  modeling  of  the  dominant 
turbulence  component  leads  again  to  dependence  as  a”3v£u)”s  (e.g. 
see  Eq.  27).  Though  it  would  be  extravagant  to  claim  robustness 
for  primitive  modeling  that  entails  so  many  uncertain  elements, 
no  definite  possible  reason  for  contrary  behavior  has  been 
uncovered. 

Some  suggestion  of  the  dependence  of  wall  pressure  at  non- 
negligible  wavenumbers  has  been  attempted  by  use  of  the  terms 
of  next  order  in  the  formal  expansion  in  terms  of  boundary-layer 
sources  (see  Eq.  32),  assuming  the  same  type  of  model  for  the 
dominant  turbulent  source  field.  A  term  proportional  to  6/a 
is  obtained  but  only  at  wavenumbers  so  low  that  the  zero-wavenumber 
part  would  predominate  there  (see  Eq.  36  with  k+<6~1  and  compare 
Eq.  27) .  A  term  is  also  obtained  that  differs  from  the  correspond¬ 
ing  term  «k2  in  the  current  model  (9a)  by  a  factor  ~  (k+a)”2 
(>>1)  at  small  k+a  (see  Eqs.  39,  7) .*  The  formal  expansion  used, 
however,  nominally  requires  that  |k|6<<l  as  well  as  DJu<x/a)«  1 ,  so 
that  this  modeling  of  the  term  of  next  order  in  k  may  be  invalid 
where  departure  from  the  zero  order  is  appreciable  at  all. 


♦The  weakly  frequency-  and  wavenumber -dependent  factor  ~|C’|2  in 
(39)  (see  Eq.  32a)  is  also  a  new  feature  that  will  persist  in 
the  corresponding  contribution  more  generally. 
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In  view  of  this  latter  limitation  it  seems  useful  for 
further  semiempirical  modeling  to  seek  the  formal  expansion  to 
the  subject  order  where  the  condition  M<$<<1  is  relaxed  (but 
still  MUoo/uj<<  1 )  .  The  corresponding  generalization  was  carried 
out  in  Ref-  3,  Eq.  35,  in  the  instance  of  planar,  incompressible 
flow.  In  further  work  the  same  may  be  attempted  for  the  axi- 
symmetric  component  in  incompressible  cylindrical  flow.* 


APPENDIX  1.  INDICATIVE  MODELING  OF  TURBULENCE  SPECTRA  IN 
CYLINDRICAL  COORDINATES 

As  suggested  in  Sec.  4,  and  underlying  modeling  pursued 
there,  model  spectra  are  constructed  here  in  cylindrical  co¬ 
ordinates  without  explicit  reference  to  a  cylindrical  wall  as 
the  physical  generator  of  the  turbulence.  Consider  first  a  cross- 
spectrum  E  (r,r 1  ,<p,k,ui)  between  amplitudes  of  a  quantity  for 
axial  (x)  wavenumber  k  and  frequency  to,  evaluated  at  position 
vectors  r,r'  in  the  plane  normal  to  x  with  polar  coordinates 
(r,0),  (r’,0+<j>).  The  random  process  assumed  is  statistically 
axisymmetric  in  the  normal  plane,  so  that  E  in  fact  depends  on  the 
angular  coordinates  only  via  their  difference  (see  sketch) . 

E  may  be  decomposed  into  angular  harmonics. 


(Al-la) 

(Al-lb) 


E  (r,r '  ,<t>,k,w)  =  I  e^^E  (r,r',k,to), 
m=-°°  ^  m 

Em(r,r '  ,k,io)  =  (2tt)  " 1 1  d^e'^E  (r,r  '  ,<{>  ,k,to)  . 

-IT 


*By  analogy  with  the  planar  result,  it  is  expected  that  the  thus 
generalized  form  of  Eq.  32  will  involve  a  factor  in  the  integrand 
that  decreases  with  ky. 
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Suppose  initially  that  the  process  is  homogeneous  and  isotropic  in 
the  normal  plane;  then  the  dependence  of  E  may  also  be  expressed  by 

E  (r,r  '  ,<j)  ,k,u>)  ■  E (R($) ,k,w) , 

R  (4> )  =  |r'-r|  =  (r2+r '  2-2rr  'cos<J>)  **  (Al-2) 

For  representative  modeling,  as  discussed  in  the  text,  suppose 
that  E  depends  on  k  and  to  only  via  k+,  defined  in  Eq.  7,  and  on 
a  single  length  scale  A  and  velocity  scale  v* .  For  definiteness, 
consider  the  cross-spectrum  of  the  axial  component  of  fluctuating 
velocity:  E  =  E^,  say.*  Then 

E  (r,r 1  ,<j>,k,(o)  =E  (R(<|>)  ,k,u»)  -  v*A2F  (k.A,R(<j>)/A)  , 
ii  ii  n+ 

(Al-3) 

where  F  is  a  dimensionless  function  of  the  given  arguments. 

By  Fourier  transformation  of  E  also  over  the  coordinates 

i  i 

of  separation  in  the  normal  plane,  a  four-dimensional  spectral 
density  may  be  defined  that  depends  on  the  three-dimensional 
wavevector  k  only  via  its  axial  and  normal  components  k  and  K,  say. 
This  spectral  density  is  then  expressed  inversely,  using  (Al-3) ,  by 

00 

E  (k,K,u)  =  C2tt) “ 2  (2tt)  f  dRRJ  (KR) E  (R,k,u>) 
ii  Jo  o  11 

-  v*A**G  (KA,k  A) , 

1  1  + 

00 

G  (KA ,k  A)  =  (2TT)"1!  dzzJ  (KAz)F  (k,A,z)  .  (Al-4) 

ii  +  J.  o  11  + 


*It  would  be  expected  that  the  spectral  density  of  one  component  or 
function  of  velocity  may  be  more  nearly  a  function  of  k+  than  another 
and  that  the  value  of  the  velocity  scale  factor  h  in  Eq.  7  for  k. 
that  renders  this  most  nearly  true  will  also  depend  on  the  identity 
of  the  component  or  function  (cf.  Ref.  21,  Sec.  2.3). 
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P  and  G  in  Eqs.  Al-3  and  Al-4  may  be  related  to  the  corre- 
1111 

sponding  space-time  correlation,  which  under  present  assumptions 
may  be  written,  for  axial  separation  x,  normal  separation  R, 
and  time  delay  r, 

W  (x,R, t)  ■  vjf (x  /A,R/A) , 

11  + 

X.  »  [  (X-U_T)  2+  (hv*T) 2  ]**  ; 

T  C 


explicitly, 

00 

F  (k.  A  ,  R/A)  *  (2n) *1h"1 fdzzJ  (k  Az) f (z,R/A) , 

ii  +  J0  •  * 

00  00 

G  (KA,k.  A)  =  (2tt) -2h” 1  f  dzzJ  (k,  Az)  f  dCCJ  (KAe)  f  (z,  C)  • 
n  +  J0  0  jo 

Likewise,  for  representative  modeling  suppose  that  in  the 
spatial  correlation  function  the  axial  and  normal  components  of 
separation  add  in  quadrature;  then  (Al-5)  assumes  the  still  more 
special  form 

W  (x,R,x)  *v2f(r/A), 

ii  + 

r.  =  [  (x-u_t)  2+(3R)  2+(hv*x) 2 ]**  ,  (3  *  const.) 

T  C 

and  the  dimensionless  four-spectrum  G^  of  (Al-4)  assumes  the  dual 
form 

00 

G  (K.  A)  -  (2ir)  "2  (hB)  -1  f  dss3  (K  As)  “lJ  (K  As)  f  (s)  , 

11  +  Jo  1 

K .  -  (ki+B^K2)15  =  [<w-u  k)2/(hv*)S^k2  +  S•-^K2],5  . 

A  convenient  class  of  such  correlations  and  spectra  that,  by 
variation  of  its  parameters  (A  ,b^  4^,  v),  encompasses  a  broad 
range  of  plausible  functional  behavior  is  defined  by  assuming  in 
ftl-7)  that 


(Al-5) 


(Al-6) 


(Al-7) 


(Al-8) 
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f<p)  -  A  {p2+b2)’v/V  (b  (p^c2)**)  .  (Al-9) 

o  1  v  o  1 

Performing  the  required  transforms  by  Ref.  22,  p.  72,  (35)  yields 
for  corresponding  spectra  (see  Eqs.  Al-3,  Al-4) 

F  (k  i,R/4)  *  (2tt)-1A  h“1b"V(62R2/A2+b2)  (1_v)/2(k2A2+b2)  (v-1)/2 
11+  0  0  1  +0 

X  K,  .  ((82R2/A2+b2)‘S(k2A2+b2)Js)  ,  (Al-10) 

V-l  1+0 

G  (K  A)  *  (2tt)  “2A  h-1g-2b"Vb2“v(K2A2+b2)”1+V/2K  ~(b  (K2A2+b2)>I) 

11+  0  Ol  +  o  V-Z  1  +  0 

(Al-11) 

It  is  recalled  that  K_v(z)  =  Kv(z);  thus  the  class  of  functions 
(Al-11)  for  the  four-spectrum  is  the  same  as  the  class  (Al-9) 
for  the  space-time  correlation.  This  is  the  same  class  of  func¬ 
tions  introduced  in  Ref.  11,  Eq.  36. 

For  b^O,  the  spectral  dependence  (Al-11)  becomes  power-law: 

G  «  (K2A2+b2)"(2”v)  . 

11  +  o 

If  also  v  *  0,  it  is  noted,  the  result  (Al-4)  for  Elt  when  K+A>>bo 
becomes  scale- independent,  i.e.  the  dependence  on  A  cancels. 

To  achieve  the  fit  to  the  experimental  results  of  Ref.  18 
discussed  in  Ref.  11,  p.  43,  referring  to  a  planar  wall,*  with 
the  present  scale  A  taken  proportional  there  to  K  =  (yy,)>s,  by 
inspection  of  Eqs.  40,  36,  and  the  choice  v  «  1  there,  the  corre¬ 
sponding  assumption  here  is  v  ■  1.  Other  parameters  of  that  fit 
imply  also  that  b  "0  and  b  A  ■  (2yy')J*  (in  the  planar  limit). 


♦The  fit  to  this  experiment,  as  noted  in  Ref.  11,  does  not  involve 
the  (weak)  temporal  dependence  in  the  convected  frame  nor  hence 
validate  the  modeling  assumption  of  space-time  isotropy  (within 
a  scale  factor) . 
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The  choice  of  a  similar  model  for  the  cylindrical  case  was  intro¬ 
duced  following  Eq.  17  of  the  text.  The  present  choice  of 
parameters  will  be  referred  to  here  as  the  "standard  model." 

The  models  considered  in  this  appendix  implicitly  assume 
the  scale  A  to  be  independent  of  the  spatial  points  involved 
in  the  correlation  referred  to.  In  the  actual  instance  of  wall 
turbulence,  as  discussed  in  the  text,  dependence  of  A  on  y,y' 
is  simply  assumed  afterward  in  the  definition  of  the  intended 
model  from  the  formally  constructed  one,  as  implied  in  the 
preceding  paragraph. 

For  the  above  standard  model,  the  cross-spectrum  (Al-3)  has 
the  form 

^  (r,r '  ,<f>  ,k,ui)  =  A'v*A2K(i  (k+(62R2+b*A2),s)  .  (Al-12) 

We  proceed  by  rough  approximations  to  simplify  the  form  of  the 
model  as  used  in  the  text  and  to  explore  its  characteristics. 

By  Eq.  Al-lb  the  m-th  harmonic  cross-spectral  component, 

E®^  (r,r '  ,k,u>) ,  is  given  by 

-IT  5 

E®  (r,r',k,oi)  -  (2ir)-1A'v*A2f  d^e”1®^  (k,  (82R2  (<|>)+b2A2)  **) 

11  i  o  +  1 

(Al-13) 

with  R($)  as  in  (Al-2) . 

For  r«r',  one  has  R($)  *  2r|sin^4>|.  For  r,r'  not  necessarily 
equal  but  $<1,  one  has  roughly 

R ( 4> )  ~  [  (r  '-r)  2+rr  '<j>2  ]**  ;  (Al-14) 

in  the  most  contrary  case  where  the  two  points  are  diametrically 
opposed  ($-ir)  ,  this  yields 

R($)  ~  (r ’+r)2  +  (ir2-4) rr '  , 
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while  the  exact  result  is  R(u)  =  (r'+r)2.  For  indicative,  order- 
of-magnitude  scaling  estimates,  as  in  the  text,  it  appears  that 
approximation  (Al-14)  will  suffice. 

Eqs.  Al-13  and  Al-14  yield  for  the  axisymmetric  harmonic  m«0 

E  °(r,r '  ,k,a»)  =  ff_1A'v*A2I  , 

11  o 

.IT 

I  -  I  d$K  (k  (A2+82rr *<|>2)  **)  , 
o  Jo  o  +  r 

A2  =  02 (r '-r) 2+b|A2. 

By  use  of  limiting  forms  of  (z)  for  small  and  large  arguments, 

I  can  be  estimated  in  pertinent  distinct  domains.  In  an 

0  L 

approximation  where  in[k+B(rr')  ]  is  not  distinguished  from  a 
constant,  the  result  may  be  written  simply  as 

1,  k+Ar<l,k+(rr ' 

E°  ^  (r,r '  ,k,u)  ~  v*A2*^  [k+(rr  • ) **63  ,  k+Ar<l,k+(rr ')  SrB^l 

i[k+(rr'),sa]'*exp(-k+Ar)  ,  k+Ar>l 

(Al-15) 

I 

This  estimate  underlies  and  exemplifies  form  (15)  in  the  text. 

Referring  for  a  moment  to  arbitrary  harmonics  m,  it  may  be 
expected  from  Eq.  Al-lb  or  Al-13  that  E^  (r,r '  ,k,u»)  will  remain 
substantial  relative  to  the  value  for  m*0  for  |m|<l/A$,  where  A<|>  is 
the  halfwidth  in  4>  over  which  E  (r,r',$,k,u)  remains  relatively 
appreciable.  Hence  the  point  spectrum  (r'»r)  is  estimated  by 


E  (r,k,oj) 
i  i 


£  E™  (r,r,k,u>) 
m  1 1 


ME  9  (r,r ,k,u>)  , 
i  l 


M 


A<|>>1 
A4><  1 . 


For  the  standard  model  (Al-12)  ,  with  r'-r,  one  has  A^-b^/Sr. 
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We  consider  now  briefly,  in  place  of  the  standard  model, 

one  where  in  (Al-9) -{A1-13J  b  -*■<)  and  u  =  l-v>0.  By  (Al-10) 

1 

E  (r,r '  ,<j),k,o))  +  Bv*A2(8R/A)y (k2A2+b2)“p/2K  ( (k2 A2+b2 ) ^BR/A)  . 

11  +0U+0 

(Al-16 ) 

In  this  case,  in  place  of  estimate  (Al-15)  one  obtains 

(  1*  kA|r'-r|@a,  k^rrD^Bsl 

r*A2  (kAA)“2yx  /  [kA (rr ' )  **  3]” 1 ,  kA|r'-r|S<l,  kA  (rr ' )  ^l 
[k  (rr  ' ) 55  g]‘ 1  [kA|r'-r|  g]yexp(-6kAlr '-r| ) , 
"  '  kA|r'-r|B>l 

(Al-17) 

where  kA  =  (k£+b2A“2)  \  The  halfwidth,  for  r'=r,  is  given  by 
A(j»~l/6kAr. 

Estimates  (Al-15)  and  (Al-17)  can  be  subsumed  in  an  expression 
valid  for  the  class  of  functions  (Al-9) - (Al-11)  more  generally. 


E  0  (r,r ' , k, to)  ~  v* 
1 1 


APPENDIX  2.  MEAN  VELOCITY  PROFILE  AND  BOUNDARY-LAYER  GROWTH  IN 
TURBULENT  FLOW  ALONG  A  CYLINDER:  COMPARISONS 
WITH  RECENT  MEASUREMENTS 

A. 2.1  Mean  velocity  profile 

Reference  was  made  in  this  memo  to  the  dependence  of  a 
mixing  length  derived  for  the  "logarithmic  range"  of  Rao's  pre- 
seciption20  for  the  mean  velocity  in  flow  along  a  cylinder.  For 
possible  future  reference  as  well,  we  therefore  append  a  compari¬ 
son  of  a  result  of  this  prescription  with  the  mean  profiles 
measured  by  Willmarth  et  al.19  on  cylinders  of  various  radii  and 
by  Leehey  and  Stellinger5  on  a  (thin)  cylinder  at  various  stream- 
wise  positions.  These  measurements  are  much  more  extensive  than 
the  earlier  ones  by  Richmond25  compared  in  this  way  previously. 2 6 
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The  Rao  prescription  is  reviewed  as  follows.  The  mean 
velocity  profile  U(x,y)  in  the  turbulent  flow  (without  pressure 
gradient)  may  be  written,  over  some  inner  range  of  wall  dis¬ 
tances  y  (limited  at  least  by  the  condition  y<6),  in  the 
functional  form 

U/v*  *  F(v*y/v,y/a)  (A2-1) 

(where  v*  may  depend  on  x) .  Rao's  hypothesis  prescribes  the 
function  F  for  the  cylinder  in  terms  of  its  limiting  form,  F(v*y/v,0), 
for  a  planar  boundary: 

F(v*y/v,y/a)  ■  F ( (v*a/v) £n(l+y/a) , 0) .  (A2-2) 

Within  the  viscous  sublayer  where  F(v*y/v,0)  ■  v*y/v,  this  pre¬ 
scription  yields 

U/v*  »  (v*a/v)S,n(l+y/a) , 

which,  apart  from  terms  of  the  order  of  (y/a)“,  correctly  repre¬ 
sents  the  result  for  a  laminar  cylindrical  boundary  layer  and 
for  the  viscous  sublayer  of  the  subject  turbulent  one.  Outside 
this  domain  (i.e.,  for  v*y/v>6  or  y/a>l) ,  prescription  (A2-2) 
apparently  lacks  any  sustainable  basis  except  as  a  plausible 
semiempirical  trial  form.27 

As  in  Ref.  26,  we  apply  (A2-2)  numerically  here  on  assumption 
that  the  profile  F(y+/0)  for  a  planar  boundary  is  adequately 
approximated  by  use  of  Squire's  transitional  profile  and  without 
regard  to  changing  dependence  in  the  outer  wake-like  region 
where  y-6: 
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r 


y+ 


0<y+<7 .4 


F(y+,0)  *  / 


(A2-3) 


|2.5£n(y+-4.9)+5.1,  y+>7.4 

Results  of  Eqs.  A2-1  to  A2-3  are  compared  in  Fig.  A2-1  with 
representative  velocity  profiles  measured  by  Willmarth  6t  <zZ.19 
on  cylinders  of  various  radii  0.01  in.  to  1  in.,  as  given  in 
their  Fig.  32.  The  assumed  profile  of  Eq.  A2-3  for  planar 
boundary  is  also  shown  along  with  that  of  Coles2®  as  given  in 
Ref.  19.  The  computed  profiles  are  seen  to  provide  a  fair 
representation  of  the  data.  More  specifically,  the  fit  is 
excellent  for  all  v*y/v  (even  such  that  y/a>>l)  at  6/a  s  16 
(with  v*a/v~200),  apart  from  apparent  inexactitude  of  the  assumed 
form  (A2-3)  in  the  transitional  region.  The  measured  dependence 
on  6/a  or  on  v*a/v  at  fixed  v*y/v,  however,  where  5<6/a<27, 
is  greater  than  that  computed. 


Computed  results  are  similarly  compared  with  profiles  mea¬ 
sured  by  Leehey  and  Stellinger5  in  Figs.  2  to  4,  as  given  in  their 
identically  numbered  figures.  The  data  in  each  figure  refer  to 
a  fixed  free-stream  speed  and  hence  roughly  equal  values  of  v*a/v 
but  to  a  wide  range  of  distances  x  along  the  cylinder  and  hence 
also  of  6/a.  The  computed  curve  shown  in  each  figure  corresponds  to 
choice  of  the  smallest  of  the  pertinent  series  of  values  of 
v*a/v  quoted  by  the  authors,  but  for  the  spread  of  values  of 
v*a/v  in  each  instance  the  range  of  variation  of  U/v*  at  fixed 
v*y/v  is  negligible— only  about  0.2  at  the  greatest  values  of 
v*y/v  and  less  at  smaller.  The  assumed  profile  (A2-3)  for  a  planar 
boundary  is  again  shown  in  Fig.  2  and  in  Figs.  2  to  4  of  the 
Coles  profile  as  given  in  Ref.  5. 


The  measured  results  display  substantial  dependence  on  6/a, 
especially  at  the  largest  vAa/v,  extending  in  this  case  even  into 
the  range  where  y/a<l.  The  calculated  results  in  each  figure  agree 
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very  well  with  the  measured  profiles  for  6/a-15  to  20. 

For  the  assuned  6-independent  planar  profile  (A2-3) ,  however, 
the  computed  profiles  do  not  depend  on  6  (at  fixed  v*a/v)  and 
do  not  capture  the  decrease  with  6/a  displayed  (at  fixed 
v*y/v)  by  the  measured  profiles.  This  dependence  on  6/a  was 
emphasized  and  discussed  under  a. different  assumption  by 
Leehey  and  Stellinger.s  The  sense  and  magnitude  of  the  dis¬ 
crepancy  at  larger  6/a  parallels  that  noted  above  in  Fig.  1 
in  the  instance  of  the  measurements  of  Ref.  19.  (See  also  Ref.  30.) 

The  Rao  hypothesis,  in  summary,  seems  to  provide  a  useful 
rough  description  of  the  mean  velocity  profile  on  cylinders 
but  does  not  encompass  the  full  dependence  on  the  characterizing 
variables  6/a  and  v*a/v.  This  limitation,  moreover,  would  not 
be  removed  by  including  the  dependence  on  the  outer  variable  y/6 
in  Eq.  A2-2  as  it  enters  the  profile  for  a  planar  boundary. 

As  a  matter  merely  of  curve-f itting,  on  the  basis  of 
Figs.  1  to  4  a  good  approximation  to  mean  profiles  outside  the 
transition  layer  is  obtainable  by  retaining  Eq.  A2-2  but  adjoining 
to  F  in  Eq.  A2-1  a  factor  M(6/a),  say,  where  M(0)  »  1  and  where 
M(6/a)  increases  sharply  to  a  broad  maximum  of  about  1.1  for 
6/a- 5  to  9,  returns  through  the  value  1  again  at  about  6/a  ■  18, 
and  decreases  to  about  0.8S  for  6/a~27  to  48.  (Cf.  Ref.  5,  Eq.  9.) 

To  retain  a  good  approximation  in  the  transition  and  sublayers, 
however,  since  the  required  correction  tends  to  unity  there,  one 
would  need  to  let  M(6/a)  depend  also  on  v*y/v  for  v*y/v<25. 

A  factor  M(6/a)  independent  of  y  outside  the  transition  layer, 
it  is  noted,  would  imply  a  corresponding  correction  factor 
l/M(6/a)  in  the  mixing  length  inferred  from  Rao's  hypothesis  (Eq.  25). 
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A. 2. 2  Boundary-layer  growth 

A  simple,  rudimentary  account  of  growth  of  the  turbulent 
boundary  layer  on  a  cylinder  can  be  given  in  the  approximation, 
closely  correct  for  the  parameter  domain  of  concern  for  towed  arrays 
where  the  wall  friction  velocity  v*  is  independent  of  distance  x. 

The  resulting  form  may  then  be  appraised  by  comparison  with 
measurements  of  6  by  Willmarth  et  at}*  and  by  Leehey  and  Stellinger, 
where  6  is  defined  by  the  condition  U(x,6)  *  0.9911^. 

First,  d6/dx  is  of  the  order  of  the  ratio  of  the  normal  (y) 

velocity  to  at  the  edge  of  the  layer,  i.e.  the  boundary  of 

the  streamlines  of  rotational  motion.  By  similarity  considerations 

(see  R.  E.  Kronauer,  Ref.  6),  this  normal  component  will  be 

nearly  proportional  to  the  rms  fluctuating  streamwise  component. 

By  measurements  of  Leehey  and  Stellinger,  discussed  by  Kronauer6 

and  referred  to  in  the  above  text,  the  latter  component,  for  wall 

distances  outside  the  transition  layer  out  nearly  to  y»6  varies 
u 

as  v*(a/r)  .  Accordingly,  we  have 

dfi/dx  =  D(v*/UjaV<a+6)*  ,  (A2-5) 

where  D  is  a  numerical  constant.  Integration  of  (A2-5)  from  a 
virtual  origin  x=0  (defined  by  extrapolation  to  6=0)  yields 

6/a  »  [ 1+ ( 3/2 ) D (v*/Oa) x/a] 2 / 3 -1  .  (A2-6) 

In  opposite  limits  this  becomes 

6/a-  [(3D/2)  (v*/Uw)x/a]2/3  ,  6/a»l 

(A2-7) 

6  -  DVa/U^  ,  6/a<<l 

(but,  contrary  to  our  assumption,  v*  is  not  so  nearly  constant  in 
the  latter,  planar  limit) . 
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A  second,  perhaps  less  convincing  approach  to  this  result 
is  based  on  considering,  as  in  the  treatment  of  drag  on  wake¬ 
generating  bodies  (e.g.,  see  Ref.  29,  Secs.  21,  36),  the  transport 
of  momenton  through  a  closed  surface  stirrounding  the  cylinder. 

This  surface  may  be  conceived  so  that  the  only  nonvanishing  con¬ 
tribution  arises  from  flux  through  a  circular  area  ir(a+6)2  of  a 
normal  plane  drawn  somewhat  aft  of  the  real  or  imagined  aft 
terminus  of  the  cylinder.  With  friction  velocity  still  approxi¬ 
mated  as  constant  over  the  length  x  of  the  cylinder  upstream  of 
this  plane,  conservation  of  momentum  implies  an  equation  for  drag: 


2iraxpv2 


-pU  •  2ir 

09 


drruav(r)  ' 


(A2-8) 


where  uav(r)  denotes  the  mean  downstream  component  of  fluid 
velocity  relative  to  the  free  stream  (thus  uav<0) .  Suppose  that 
just  downstream  of  the  terminus  this  incremental  velocity  profile 
due  to  the  dragging  of  fluid  by  the  wall  of  the  cylinder  upstream 
is  nearly  proportional  to  that  of  the  rms  fluctuating  velocity 
components  at  the  terminus,  i.e.  uav»-Cv* (a/r) ^  for  a<r<a+6 , 
where  C  is  a  constant.  Then  Eq.  (A2-8)  becomes 


a+o 

2iraxpv£  «  lirpv^O^C J  drrta/r)1*  ; 


integration  yields  again  just  the  result  (A2-6)  provided  that 
OD*1. 

To  check  (A2-6)  or  (A2-7)  against  measured  results  and  infer 
a  rough  value  for  D,  one  may  form  (in  the  instance  of  Eq.  (A2-7) 
for  6/a>>l) the  ratio  (6/a)/[ (v*/Uj (x/a) ]2/s  for  the  various  data 
sets  of  Refs.  5  and  19  and  see  if  it  is  nearly  constant  (or  at 
least  has  no  substantial  systematic  variation) .  For  the  14  data 
sets  of  Ref.  19,  Table  1,  excluding  two  (3  and  12)  rejected  as 
out  of  line  with  the  remaining  ones,  the  inferred  value  of 
the  cited  ratio  varies  by  less  than  a  factor  two  (not  monotonically 
with  S/a)  while  the  range  of  hose  diameters  encompasses  a  factor 
100.  The  mean  value  inferred  for  D  is  about  0.45.  Similarly,  for 
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the  23  possible  data  sets  of  Ref.  5,  Table  1,  accepting  the 

value  of  v*  therein  obtained  by  fitting  velocity  profiles  in 

the  buffer  layer,  the  cited  ratio  varies  by  a  factor  of  about 

two  while  the  range  of  positions  x  is  10.5.  The  mean  value 

inferred  for  D  in  this  case  is  about  0.8.  The  corresponding 

value  C  *  1.25  yields  a  result  that  u  „  =  -u  (at  wall  distances 
-k  av 

where  u«r  ^) ,  where  u  is  the  measured  rms  streamwise  velocity 
component,  as  seen  in  Fig.  2.3-6  of  Ref.  6  (based  on  Fig.  3  of 
Ref.  5) . 
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CAPTIONS  FOR  FIGURES 

1.  Root-mean-squared  longitudinal  component  of  turbulent 
velocity  vs  wall  distance  for  thick  layer  on  cylinder 
(from  Ref.  6,  Fig.  2.3-6,  based  on  Ref.  5,  Fig.  8). 

A2-1.  Measured  and  computed  mean  velocity  profiles  for  cylinders 
of  different  diameters.  Measured:  Fig.  32  of  Ref.  19; 
computed:  Eqs.  A2-1  to  A2-3  with  av*/v  from  Ref.  19, 

Table  1,  tabulated  below  (and  planar  limit  av*/v*») ; 
also.  Coles’  profile  for  planar  flow  from  Ref.  19,  Fig.  32. 


Symbol 

6/a 

av*/v 

9 

1.8 

2710 

A 

4.1 

1360 

A 

4.7 

751 

■ 

9.4 

389 

□ 

16.0 

198 

• 

27.0 

83.3 

O 

37.5 

46.1 

Measured  and  computed 

mean  velocity  profiles  for  a  cylinder 

at  different  longitudinal 

positions 

.  Measured:  Fig.  2  of 

Ref.  5  (see  table  below); 

computed: 

Eqs.  A2-1  to  A2-3 

for  av*/v  »  57.9  (and 

for  planar  limit  av^/v-^) ;  also. 

Cole' a  profile  for  planar 

flow  from 

Ref.  5,  Fig.  2). 

Symbol 

6/a 

a  v*/v 

4 

25.5 

57.9 

A 

26.4 

59.1 

■ 

25.6 

59.9 

□ 

28.1 

58.7 

• 

29.2 

59.8 

O 

41.4 

61.3 
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CAPTIONS  FOR  FIGURES  (Cont'd.) 

A2-3.  Profiles  as  in  Fig.  A2-2  but  from  Fig.  3  of  Ref.  5  for 
parameters  tabulated  below, and  computed  profile  for 
av*/v  *  81.9 


Symbol 

X  € 

a 

e 

A 

A 

a 

a  o 

6/a 

15.3  20.0 

17.3 

24.5 

26.8 

29.0 

31.4 

44.1  38.4 

a  v*/v 

83.3  81.9 

83.4 

83.4 

84.1 

85.7 

86.9 

88.1  90.6 

A2-4  Profiles  as  in  Fig 

.  A2-2 

but  from  Fig 

.  4  of 

Ref. 

5  for 

parameters  tabulated  below, and 

computed  profile  for  av*/v  = 

Symbol 

<d  a 

e 

A 

A 

a 

□ 

O 

6/a 

18.3  16.3 

23.3 

27.4 

26.0 

24.5 

45.4 

47.5 

av*/v 

103.2  102.8 

105.3 

108.1 

116.2 

117.3 

115.0 

115.0 
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